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Abstract — In this paper a novel filtering design intended for 
the impulsive noise removal in color images is presented. 
The described scheme utilizes the rank weighted cumulated 
distances between the pixels belonging to the local filtering 
window. The impulse detection scheme is based on the 
difference between the aggregated weighted distances assigned 
to the central pixel of the window and the minimum value, 
which corresponds to the rank weighted vector median. If the 
difference exceeds an adaptively determined threshold value, 
then the processed pixel is replaced by the mean of the 
neighboring pixels, which were found to be not corrupted, 
otherwise it is retained. The important feature of the described 
filtering framework is its ability to effectively suppress 
impulsive noise, while preserving fine image details. The 
comparison with the state-of-the-art denoising schemes 
revealed that the proposed filter yields better restoration 
results in terms of objective restoration quality measures. 

Index Terms — color image processing, noise reduction, image 
enhancement 

I. Introduction 

Quite often the quality of color images is degraded by 
various types of noise, whose removal is required to enable 
the success of further steps in the image processing pipeline. 
Noise, arising from a variety of sources, is inherent to 
electronic image sensors and therefore the noisy signal has 
to be processed by a suitable filtering algorithm that reduces 
the noise component, while preserving original image details 
[1-3]. 

In this work we focus on a special kind of distortion, 
which is introduced by impulsive noise caused by 
malfunctioning sensors, faulty memory locations in hardware, 
aging of the storage material or transmission errors due to 
natural or man-made processes [4-6] . 

The color image will be treated as a two-dimensional ma- 
trix consisting of pixels y = fx „ y „ r 1, where the 
index u = 1,...,N indicates the pixel position on the image 
domain and v denotes the number of channels. In the case 
of standard color images v=3. The vector components x ., 
for k = 1,2,3 represent the RGB color channels values. 

Generally, filtering operators work on the assumption that 
the local image features can be extracted from a small image 
region centered at pixel x , called a sliding filtering window 
and denoted as w. Thus, the output of the filtering operation 
will depend only on the n samples contained within the 
filtering window centered at x , which will be also denoted 
for convenience as , so that W — {x^,x?, ...,x n }. 

The majority of the nonlinear filters used for impulsive 

©2013 ACEEE 44 
DOI:01.IJSIR4.3.1255 



noise removal in color images is based on the order statistics 
[7] . These filters perform the vector ordering of the set of 
pixels from the filtering window to determine the filtering 
output. The widely used reduced vector ordering is based 
on assigning a dissimilarity measure to each color pixel from 
the filtering window [1,8,9]. The aggregated dissimilarity 
measure assigned to pixel x is defined as 



Ri = P ( X i> X j) • X i> x j E 



W, 



(1) 



/J 



where p(') denotes the distance between pixels in a given 
color space. The values of R ., (i = l,...,n) are then sorted and 
the vectors x p x 2 , x n are correspondingly ordered 

< R f7l) =? x <X) < < x {n> (2) 



where < denotes the order relation between vectors and 
R , , denotes the k-th smallest value of R. 

(K) 

Many denoising techniques define the vector x (1 . in (2) 
as their output, since vectors that diverge significantly from 
the samples of w appear in the higher indexed locations in 
their ordered sequence. 

One of the most widely used noise reduction techniques 
is the Vector Median Filter (VMF), whose output is the vector 
x.j. from W, for which the sum of distances to all other vectors 
is minimized 



X(±y = argmin sew £/=i|f - Xj\l (3) 
where 1 1 ■ 1 1 denotes Euclidean norm. The VMF is the most 
popular operator intended for smoothing out the impulses 
injected into the color image by the noise process. This filter 
is very efficient at reducing the impulses, preserves sharp 
edges and linear trends. However, the drawback of the VMF 
and other filters based on vector ordering lies in introducing 
excessive smoothing, which results in an extensive blurring 
of the output images. Thus, the filters based on vector order- 
ing do not preserve fine image structures, which are treated 
as noise and therefore generally show a tendency to the 
blurring of image details and generation of color artifacts. 

In order to alleviate the problem of image smoothing 
various switching filters, that replace only the corrupted 
pixels have been proposed [10-13]. The efficiency of a 
switching filter depends both on the quality of the impulse 
detection scheme and on the applied restoration framework, 
which replaces the detected impulses with estimates derived 
from the samples belonging to a local processing window. 

In this paper a novel switching filter is proposed. The 
main advantage of the proposed approach is its ability to 
suppress the noise component, while preserving fine image 
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details. The structure of the filter is based on the reduced 
ordering statistics and is characterized by a low 
computational complexity, which enables the adoption of the 
new technique in real-time applications. 

II. Rank Weighted Filtering Design 

The reduced ordering schemes are based on the sum of 
the dissimilarity measures between a given pixel and the 
samples from the local filtering window w. In this way, the 
output of the vector median filter is the pixel whose average 
distance to other pixels is minimized. 

The distances Pij = P{%i, ^;)between the pixel xand all 
other pixels w belonging to can be arranged into a sequence 

Pi±,Pi1* - -.Pin -» P;(l> Pi(2> ■■■.Pi(ri), (4) 

and the ranks of the ordered distances can be used for the 
evaluation of the aggregated distances in (1). 
If stands for the rank of a given distance, then Pi(r) vill 
denote the corresponding distance value and instead of the 
aggregated distances in (1) a weighted sum of distances, 
utilizing the distance ranks, can be composed 



(5) 



where /f r) is a decreasing weighting function of the distance 
rank r. Then, the rank weighted sums of distances D can be 
sorted and a new sequence of vectors is obtained 

B(1},B(2}, -•■ » ^%i>%2>--*%n> (6) 

where the vector ^(i) is the output of the Rank Weighted 
Vector Median Filter (RWVMF). Applying a step-like func- 
tion 

jr/ys _ fl. f° r r £a, ex < n, (J) 
' lO, otherwise, 

the Sharpening Vector Median Filter (SVMF) presented in 
[14] is obtained. 

Extensive experiments revealed that satisfactory denoising 
results are achieved using monotonously decreasing function 
f(r) = 1/r [15-17]. The weighting function decreases the 
influence of large distances introduced by the outliers injected 
by the noise process, which enables to efficiently remove 
the impulsive noise while enhancing the image edges. 



proposed method. It can be observed that the map of the 
detected noise correlates very well with the real contamina- 
tion measured as the Euclidean distance between the origi- 
nal and corrupted pixels in the RGB color space normalized to 
the range < 0,255>. Figure 3 exhibits the correlation be- 
tween the real and the detected impulsiveness of the image 
pixels using the test images Goldhill and Parrots. The re- 
sults reveal a high correlation between the real and detected 
noise. 




c) map of the real noise d) map of Dj — D m 

Figure 1. Illustration of the noise detection efficiency using the 
Parrots image corrupted by impulsive noise with p=0.1 intensity. 
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III. Adaptive Switching Filter 

In order to decide whether a pixel of a color image is cor- 
rupted by impulsive noise, the difference between the cumu- 
lated weighted distance D / assigned to the central pixel of 
the filtering window and the value of D corresponding to 
the rank weighted vector median filter output can be used 
[15-17] and the strength of the impulsive contamination can 
be estimated as the difference between Z^and D , 

Figures 1 and 2 show examples of the detected noise us- 
ing parts of the color test images Parrots and Goldhill cor- 
rupted by impulsive noise, whose intensity p denotes the 
percentage of corrupted pixels. The visual comparison of the 
real and detected noise confirm the good efficiency of the 
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c) map of the real noise i) map of 1 - D (T; 

Figure 2. Illustration of the noise detection efficiency using the 
Goldhill image corrupted by noise with p=0.1 intensity. 

Figure 4 shows the histograms of the corresponding values 
of D p D.j and the measure of impulsiveness D ; - D for 
different levels of noise intensity using the Parrots test image. 
As can be seen, with increasing level of noise intensity, the 
histograms are shifted towards higher values, which is 
confirmed by the mean value of D t - D ; marked with a red 
line on the plots. 

The structure of the proposed switching filter is quite 
simple. If the difference Z) ; - D , which measures the pixel 
corruption, exceeds a given threshold value , then a pixel is 
declared as corrupted by a noise process, otherwise it is 
treated as not disturbed 



if D ± — D 
otherwise, 



w> ( 8 ) 
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a) GOLDHILL, p = 0.1, p = 0.97 b) PARROTS, p = 0.1, p = 0.98 




OGoi.DHi: :./>- n.2,/> - n.'w di parrots, p o.2.p 0.99 






measure D J - D computed for all image pixels and denoted 
as pL. , (see also Fig. 4). The dependence between the average 
impulsiveness measure fj. and noise intensity level p is 
shown in Fig. 7. 
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Figure 3. Illustration of the correlation between the real noise and 
the output of the proposed noise detector for different levels of 
noise intensity. The correlation coefficient p is provided for each 
plot. 

where y is the switching filter output, x is the central pixel of 
the filtering window and X is the Arithmetic Mean Filter 
output computed using only the pixels found by the detector 
to be not corrupted by the noise process. If all neighbors of 
the central pixel x / of the filtering window are declared as 
corrupted, then the VMF applied to all pixels from w is taken 
as the filter output. 

Of course, the efficiency of the switching scheme is 
dependent mainly on the value of the thresholding parameter. 
If the threshold T is too low, the filter will be replacing 
uncorrupted pixels and much of the image details will be lost. 
On the other hand, if the value of T is too high, many 
corrupted pixels will pass the filter without being processed. 

As could be expected, the optimal setting of T depends 
on the contamination intensity. As can be observed in Fig. 6, 
which shows the dependence between the optimal threshold 
X /alue and the noise contamination level p evaluated using 
a set of test images shown in Fig. 5, the threshold yielding 
the best PSNR value is decreasing with increasing noise 
contamination level and does not depend significantly on 
the image structure. 

Therefore, the thresholding parameter T needs to be 
adjusted to the noise intensity level. The experiments 
performed using the set of images shown in Fig. 5, indicate 
an approximately linear dependence between the 
contamination ratio p and the mean value of the impulsiveness 

©2013 ACEEE 46 
DOI:01.IJSIR4.3.1255 



fvj 9001 
|| 4-0(1 

sa- 
il 1 







































20 -lo Go 80 UKi 13) Mu 160 l«0 2oo 220 340 



Mill 

BOO 
Kill 

■} 10 

II 































V 
















I, 












































20 III (i 


) 8 


) loo 120 1 


101) ISO 3IKI 220 240 



rl 

Q 




CI) 80 100 13) 110 160 ISO Mi 320 310 



Figure 4. Histograms of the values of Z) ; , D and D l - Z) created 
using the PARROTS image contaminated by the impulsive noise 
with different intensities. The red line marks the mean value of a 
histogram. 
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image 1 image 2 image 3 image 4 image 5 

fti' 




image 6 image 7 image 8 image 9 image 10 

Figure 5. Color test images used for the construction 
of an adaptive filter. 
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Figure 6. Dependence of the optimal threshold T on the 
contamination level p for the images shown in Fig. 5 corrupted by 
impulsive noise. 
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Figure 1. Dependence between the optimal threshold T and the 
output of the noise detector for the images shown in Fig. 5. 

IV. Comparison with Existing Techniques 

For the evaluation of the effectiveness of the proposed 
filtering technique the test images were corrupted by a 
impulsive noise modeled in such a way, that the noisy pixel 
x t = {x iL ,x i2 ,x i3 } is defined as: 

(10) 



with probability p 
with probability 1 - 



Figure 7. Dependence of the measure of impulsiveness fJL on the 
contamination level for color images depicted in Fig. 5. 

Combining these plots, we obtain the scatter plot depicted 
in Fig. 8, which reveals a roughly linear dependence between 
the optimal value of the adaptive threshold % ar, d tne value 
of pi 

t = -2.87 ■ (i + 66.5, (9) 
which can be used for the adaptive setting of the appro- 
priate threshold value enabling high efficiency of the pro- 
posed filter for low and moderate noise contamination levels. 
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1 (o t , witfi probability 1 — p 

wherein = {v iv V i2 ,V i3 }and v ik E (0,255). k - 1,2,3 . In 
this impulsive noise model, the affected pixels have corrupted 
all three channels, which take on random values from the 
interval <0,255>. 

For the measurement of the restoration quality, the commonly 
used Root Mean Squared Error (RMSE) expressed through 
the Peak Signal to Noise Ratio (PSNR) was used, as the RMSE 
is a good measure of the efficiency of impulsive noise 
suppression and correlates well with other commonly used 
restoration quality measures [1]. 

The effectiveness of the proposed adaptive technique 
was compared with a set of switching filters intended for the 
suppression of impulsive noise in color images. The following 
filters were chosen for the comparison: 

• Peer Group Filter (Fj), [18], 

• Adaptive Center Weighted VMF(F 2 ), [19], 

• Adaptive Center Weighted Directional-Distance Filter (F ), 

[20], 

• Rank-Ordered VMF(F 4 ), [21], 

• Sigma Directional Distance Filter (F^), [22]. 

Their parameters were set according to the 
recommendations provided in the appropriate references. 

The four images depicted in Fig. 9 were contaminated by 
the impulsive noise with intensities ranging from 0.05 to 0.3. 
It is worth noticing, that these images were not included in 
the set of images used to establish the dependence expressed 
by Eq. (9). The results are summarized in Tab. 1 . 

As can be observed the proposed filtering method yields 
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results significantly superior to those obtained using the 
state-of-the-art denoising methods. 




GOLDHBLL LENA PARROTS PEPPERS 

Figure 9. Test images used for the evaluation of the restoration 
results. 



Table I. Comparison of the PSNR Values Obtained When Restoring the 
Color Test Images Contaminated with the Impulsive Noise Using the 
Proposed Technique and Some Competitive Filters 



Image 


71 
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F 2 


F 3 


F 4 


F 5 


Go Id hi 11 


0. 1 


37.14 


36. 19 


35.47 


3 3.94 


3 5.54 


34.94 


0.2 


33.85 


3 2.66 


3 1.07 


3 0.53 


3 2.26 


29.57 


0.3 


31.33 


29. 14 


26.29 


26.84 


28.52 


24. 13 


Lena 


0. 1 


39.81 


38.56 


37.87 


3 6.34 


38. 12 


3 6. 07 


0.2 


36.35 


34.54 


3 1.93 


3 1.69 


33.98 


29.46 


0.3 


33.35 


30.29 


26.62 


27.03 


29.65 


23.75 


Parrots 


0. 1 


38.60 


37.32 


3 6.73 
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3 6.29 


3 6.70 
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35.55 


3 3.98 


3 1.37 


3 1.83 


3 3.09 


29.82 
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29.43 
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0. 1 


41.62 
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36.71 


35.09 


3 7.07 


34.36 
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38.45 


33.42 


30.79 


30.27 


3 2.75 


27.59 


0.3 


34.31 


28.58 


25. 10 


24.89 


27.80 


21.98 



The results summarized in Tab. 1 are confirmed by the 
subjective analysis of the results depicted in Fig. 10, which 
shows the restoration quality achieved using the proposed 
filter as compared with other reference filters. 



IV. Conclusions 

In this paper a new switching filtering design has been 
proposed. The filter is based on the weighted cumulative 
distances between pixels, which are used for the detection of 
samples corrupted by impulsive noise process. The experi- 
ments performed on test images indicate a high efficiency of 
the proposed filtering design. Moreover, the new filter has a 
low computational complexity and simple structure, which 
makes it attractive for real-time applications. 
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